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Abstract 

The existence of extra dimensions allows the possibility that the fundamental scale 
of gravity is at the TeV. If that is the case, gravity could dominate the interactions of 
ultra-high energy cosmic rays. In particular, the production of microscopic black holes 
by cosmogenic neutrinos has been estimated in a number of papers. We consider here 
gravity-mediated interactions at larger distances, where they can be calculated in the 
eikonal approximation. We show that for the expected flux of cosmogenic neutrinos 
these elastic processes give a stronger signal than black hole production in neutrino 
telescopes. Taking the bounds on the higher dimensional Planck mass Mjj (D = 4 + n) 
from current air shower experiments, for n = 2 (6) elastic collisions could produce up 
to 118 (34) events per year at IceCube. On the other hand, the absence of any signal 
would imply a bound of Mo ^ 5 TeV. 



Introduction. We observe extensive air showers produced when a cosmic ray from 
outer space hits a nucleon in the upper atmosphere. The observed events have energies 
of up to 10 11 GeV, and their profile and distribution are consistent with a primary proton of 
extragalactic origin. In their way to the Earth these protons would interact with the CMB 
photons and produce pions: 

V + 72.7K -> A+ -> n + 7T+ (p + 7T°) . (1) 

The flux of cosmogenic neutrinos is created in the decay of the charged pions, and it will 
appear correlated with observable fluxes of nucleons and photons. 

Cosmogenic neutrinos are of great interest as probes of new TeV physics because of 
two generic reasons. First, they provide large center of mass energies. Second, the relative 
effect of new physics on the weakly interacting neutrinos is larger than on quarks or charged 
leptons, making it easier to see deviations. The signals of new physics could be detected in 
deeply penetrating air showers and neutrino telescopes. 

In particular, in models with extra dimensions and the fundamental Planck scale at the 
TeV the gravitational interactions are unsuppressed in the transplanckian regime. The 
possibility of black hole (BH) formation |3j by cosmogenic neutrinos has been discussed in 
several papers 0]-jH]- Here we will study the gravitational interaction at larger distances, 
where it can be calculated using the eikonal approximation QUI . This approximation 
involves linearized gravity and is not affected by the uncertainties in the cross section for 
BH formation. After discussing the bounds on the gravitational scale from air shower ex- 
periments, we will show that these elastic processes are more frequent than BH formation 
in neutrino telescopes. At large impact parameters the neutrino interacts, loses a small 
fraction of energy and keeps going. Telescopes could detect these processes because they 
are sensitive to events of energy three or four orders of magnitude below the typical energy 
(around 10 9 GeV) of the cosmogenic neutrinos. 

TeV gravity. The simplest picture of TeV gravity includes only two free parameters: the 
value of the higher dimensional Planck scale Mjj, and the number n of compact dimensions 
where gravity propagates. A third parameter, the (common) length 2ttR of the n dimensions, 
could be deduced from the 4d Newton constant: 

G D = (2kK)»G n = . (2) 

At processes below the model-independent signature of extra dimensions is graviton 
emission. The amount of energy radiated would be proportional to the accessible phase space 
or, in the Kaluza-Klein (KK) picture, to the number of KK modes of mass below the center of 
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mass energy. In this type of experiments for a given n one sets bounds on R and then deduces 
the limits on Md- From collider experiments one obtains Md > 1.4 (1.0) TeV forn = 2 (> 3) 
[TT] . whereas from SN1987A the bounds go up to 22 TeV for n = 2 One should keep in 
mind, however, that the gravitons emitted in the supernova explosion have a KK mass below 
»s 50 MeV. The simple picture could be modified above this energy, for example, with four 
more dimensions at R' ~ (100 GeV) -1 , which would bring the fundamental scale of gravity 
down to 1 TeV without affecting the physics in the supernova. 

The bounds obtained from transplanckian collisions are complementary in the sense that 
given n they are a direct probe of Md, and R is then adjusted in order to reproduce Gn- 
At energies above Md and impact parameters smaller than R the collision is a pure higher- 
dimensional process independent of the compactification details that fix the value of the 
effective Newton constant. The transplanckian collision does not see that the extra dimen- 
sions are compact, they could be taken infinite with no effect on the cross section. 

Neutrino-nucleon cross section. The TeV gravity model should be embedded in a 
string theory, which would relate Md with the string scale Ms- In the simplest set-up |T3] 
the standard model (SM) fields (open strings) would be attached to a 4d brane, whereas 
gravity (closed strings) would propagate in the whole Dd space. In this case 

M n D +2 = ^M£+ 2 , (3) 

with g the string coupling. The transplanckian regime corresponds then to energies above 
the string scale, where any tree-level amplitude becomes very soft. In the ultraviolet string 
amplitudes go to zero exponentially at fixed angle and, basically, only the forward (long 
distance) contribution of the graviton survives. This is precisely the regime where the eikonal 
approximation is valid. 

Let us consider the elastic collision of a neutrino and a parton that exchange .D- dimensional 
gravitons (see j3 [TU] for details). The eikonal amplitude A e ik(s, t) resums the infinite set 
of ladder and cross-ladder diagrams. It is reliable as far as the momentum carried by the 
gravitons is smaller than the center of mass energy or, in terms of the fraction of energy 
y = (E v — E' u ) / E u lost by the incoming neutrino, if y — — t/s <C 1 (s and t refer to the Man- 
delstam parameters at the parton level). In this limit the amplitude is independent of the 
spin of the colliding particles. Essentially, A e ik is the exponentiation of the Born amplitude 
in impact parameter space: 

Ae ik (s,t) = jjd 2 b e iqb (e^ - l) , (4) 

where x( s , b) is the eikonal phase and b spans the (bidimensional) impact parameter space. 
x(s, b) can be deduced from the Fourier transform to impact parameter space of AB orn (s, t). 
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Our Born amplitude comes from the t-channel exchange of a higher dimensional graviton: 



s 2 f d" qr 



A Born - M „ +2 J t _ ^ , (5) 

where the integral over momentum qx along the extra dimensions (equivalent to the sum over 
KK modes) gives an UV divergence. The magic of the eikonal amplitude is that it will be 
well defined despite we obtain it from an UV dependent Born amplitude: the contributions 
from large qr introduce corrections to the phase x{ s > b) only at small b (~ 1/qr), but this 
small b region gives a negligible contribution to A e ik in the transplanckian regime. 

From A e ik one obtains x(s, b) = (b c /b) n , with 

w = Ml!r W -J- (6) 

2 \2J M 2 D +n ' K ) 

The amplitude in Eq. (J3J) can OCni then be written as A e ik{s,q) = Amsb\ F n {b c q) , where 

F n {y) = -i r dx x J (xy) (e lx ' n - l) , (7) 



q = y/—i, and the integration variable is x = b/b c . For q < b~ x this integral is dominated by 
impact parameters around b c , and for q > b^ 1 by a saddle point at b s . As q (or y = q 2 / s) 
grows nonlinear corrections (H diagrams) become important ^Oj- For — t/s ~ 1 b s approaches 
j2] the Schwarzschild radius Rs of the system: 

' 2 v^r(^)\^ / s 



2(n+l) 



Rs= I n + 2 J \M 2 D n+4 

At b < Rs one expects an inelastic collision, with a significant emission of gravitons, and 
black hole (BH) formation. The latter possibility has been considered in several analyses 
where is it also shown that a number of factors (angular momentum, charge, geometry 
of the trapped surface, radiation before the collapse) make a quantitative estimate difficult. 
In particular, the higher curvature corrections discussed in [TH] could affect the evolution of 
the collision after the horizon has formed, making the simple picture of single BH production 
and subsequent Hawking evaporation unlikely. 

The differential cross section da e ik/dy grows as y decreases j^J- For example, taking 
M D = 1 TeV and E u = 10 10 GeV, for n = 2 (6) it is a factor of 265 (62) larger at y = 10" 3 
than at y = 0.1. The small y region corresponds to long distance processes where the 
neutrino interacts with a parton and transfers only a small fraction of its energy. This 
region is less important for a larger number of extra dimensions, since then gravity dilutes 
faster and becomes weaker at long distances. On the other hand, values of y close to 1 mean 



4 



shorter distance interactions. For example, we obtain that for y = 0.5 a 52% of the vN 
eikonal cross section comes from impact parameters b < Rs for n = 2 (or a 71% for n = 6). 
We will then separate two types of transplanckian (s > M D at the parton level) processes: 

(i) Inelastic processes where the neutrino interacts with a parton at distances b < Rs- 
The cross section for these processes, obh = TtR s , would include BH formation and hard 
scatterings with important graviton emission where the neutrino loses most of its initial 
energy. 

(ii) Elastic processes where the neutrino transfers to the parton a small fraction y of its 
energy (we take y ma x = 0.2), and keeps going. We use the eikonal approximation to describe 
these processes. They are dominated by impact parameter distances larger than Rs and 
thus non-linear effects and graviton emission are expected to be small. 

Air showers from cosmogenic neutrinos. The flux of cosmogenic neutrinos depends 
on the production rate of primary nucleons. It will appear correlated with proton and photon 
fluxes that should be consistent, respectively, with the number of ultrahigh energy events at 
AGASA and HiRes [T] and with the diffuse 7-ray background measured by EGRET 1.6 . We 
will consider two neutrino fluxes described in The first one saturates the observations 
by EGRET, whereas for the second one the correlated flux of 7-rays contribute only a 20% 
to the data, with the nucleon flux normalized in both cases to AGASA/HiRes. The higher 
flux predicts 820 down-going neutrinos of each flavor with energy between 10 8 GeV and 
10 11 GeV per year and km 2 , versus 370 for the lower one. 

AGASA and Fly's Eye are able to detect efficiently penetrating air showers of energies 
above ~ 10 10 GeV @]. When a cosmogenic neutrino of energy E u <: 10 10 GeV enters the 
atmosphere it can experience the two types of processes described above. Short distance 
collisions may produce a BH, whose thermal evaporation would start an air shower of energy 
(for negligible graviton emission) up to 0.8E U . On the other hand, if the neutrino suffers a 
long distance collision it will transfer to the parton a fraction y of energy up to y ma x = 0.2. 
These processes, well described by the eikonal approximation, start hadronic showers of 
energy yE v . 

We have calculated the combined number of events at AGASA and Fly's Eye as a function 
of Mjy (we take the exposures to penetrating showers from jlj). Within the SM one expects 
0.03 (0.009 for the lower flux) deeply penetrating showers started by a neutrino, with a 
background of 1.72 events from hadronic showers. In these experiments 1 event passes all 
the cuts, which implies |3] an upper bound of 3.5 neutrino events at 95% CL. For the higher 
flux and n = 2 we obtain 3.5 events (2.1 BH and 1.4 elastic) if Md = 1.0 TeV, whereas for 
n = 6 we have 2.6 BH plus 0.9 eikonal events if Md = 1-5 TeV. Neglecting the short distance 



5 



collisions, for n — 2 (6) we obtain 3.5 elastic events for Mb = 0.7 (0.9) TeV. For the lower 
flux, the 3.5 events limit would be obtained (from elastic processes only) for Mo = 0.4 TeV 
if n = 2 or for M D = 0.5 TeV if n = 6. 

TeV gravity events at IceCube. IceCube ^H] is a large scale (km 3 ) neutrino telescope 
currently under construction in the Antarctic ice. Its center is at a depth of 1.8 km, which 
implies that if <j u n < 0.01 mb neutrinos can reach it vertically with no previous interactions, 
whereas if o v ^ < 0.0001 mb they could also reach it horizontally after crossing 150 km of 
ice. The detector is sensitive to hadronic showers of energy E s h > 500 TeV. 

To be definite, let us consider a cosmogenic neutrino of energy E u = 10 10 GeV for 
Mr, = 1 TeV and n = 2 (6). The probability that the neutrino survives to reach the 
detector from a zenith angle 9 Z is 

P surv = exp[-X(e z )aN A }, (9) 

where X{6 Z ) ~ pi ce L{6 z ) is the column density of material (L(6 Z ) is the length of the column 
in ice) in its way to the detector, and a = (Jbh + ^sm is the inelastic cross section (we do not 
include a e ik because the elastic processes with y < 0.2 introduce a negligible distortion in the 
energy of the neutrinos that reach the detector). For cos8 z ~ 0.11 (0.45) the length L(0 Z ) is 
equal to its mean free path L = 1/ (pi ce aNA_), therefore the neutrino should typically reach the 
detector from smaller angles. Within the SM this angle would go up to cos 9 Z « —0.03. Once 
in the detector, the probability that the neutrino experiences a short distance interaction is 
given by 

Pi? = 1 - exp[-L p ice a BH N A ] , (10) 

where L ~ 1 km is the linear dimension of the detector. We obtain values of P^f 1 from 
0.06 for n = 2 to 0.2 for n = 6. On the other hand, the probability P^ that the neutrino 
interacts elastically, loses a fraction of energy between y min = (500 TeV)/E u = 5 x 10~ 5 
and y m ax = 0.2, and starts an observable hadronic shower, can be read from the expression 
above just by changing <j B h by da eik /dy integrated between 5 x 10~ 5 and 0.2. We obtain 
that P^t goes from 0.4 for n = 2 to 0.3 for n = 6. These probabilities only change a 1% if 
we take ymaxl^ or 2y max . Within the SM the probability that the neutrino starts a shower 
inside the detector is just P^f = 0.002. 

Given a cosmogenic neutrino flux $„, the number N s h of shower events at IceCube can 
be estimated as 

N sh = 2-kAT J d cos 9 Z J &E u -^-P surv P int , (11) 

where the sum goes over the three neutrino and antineutrino species, A w 1 km 2 is the 
detector's cross sectional area with respect to the v flux, and T is the observation time. For 
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Figure 1: Number of (eikonal and SM) shower events per year at IceCube for the two 
cosmogenic neutrino fluxes. 



the higher (lower) flux, within the SM we obtain 1.4 (0.5) hadronic or electromagnetic events 
(muons and taus do not shower) per year above 500 TeV. 

We now consider values of Mp above the bounds obtained from the absence of penetrating 
air showers and calculate the number of events per year at IceCube. Taking the higher flux, 
for n = 2 (6) we obtain a maximum of 118 (34) elastic events versus just 20 (24) short 
distance events. The elastic events correspond to soft processes with y < 0.2. In particular, 
for n = 2, 95 of the 118 showers have an energy below 10 s GeV, versus just 0.08 of the 20 
showers from short distance interactions. For n = 6, 18 of the 34 eikonal events are below 
10 8 GeV. 

In Fig. we plot the number of elastic events per year as a function of Mo for the 
two cosmogenic fluxes. We find that IceCube could detect TeV gravity effects above the SM 
background for Mjj up to approximately 5 TeV. The characteristic signature would be always 
a hadronic shower, as charged leptons are never produced in the gravitational interaction 
starting the shower. 

Discussion. Cosmogenic neutrinos interact with the terrestrial nucleons at center of 
mass energies ^2m^E l/ rs 10 5 GeV, so they could be used as probes of new TeV physics. 
In particular, the possibility of BH formation in models with extra dimensions has been 
entertained by several groups. These analyses are based on a geometric cross section that 
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assumes single BH production whenever the neutrino and the parton interact at impact 
parameters smaller than Rs- 

The problem with this estimate is that, despite the large energy of cosmogenic neutrinos, 
the uN cross section is dominated by the small x region, and most of the BHs produced in 
penetrating air showers or neutrino telescopes would be very light, with masses just above 
Mp. These light BHs would be very sensitive to effects like graviton emission during the 
collapse, higher curvature corrections, or non-thermal effects in the evaporation, which add 
uncertainty to the estimate. 

In this paper we have analyzed a different type of signal. It is produced when the neu- 
trino interacts elastically with the parton at typical distances larger than Rs and transfers 
a small fraction y of its energy. The process is properly described by the eikonal approxi- 
mation, and its distinct experimental signature would be a hadronic shower of energy yE v . 
Electromagnetic showers would never be produced in the initial uN interaction. 

We have computed the number of penetrating air showers from soft (y < 0.2) elastic 
processes and from short distance {b < Rs) processes, and have found that the expected 
number of events at AGASA and Fly's Eye from these two types of processes is similar. We 
obtain bounds on M D that are below the ones from SNf987A for n — 2, but are similar to 
the limits from collider experiments for any value of n. 

In contrast, we have shown that these elastic interactions provide a clear and model- 
independent signal of TeV gravity that would dominate over BH production in neutrino 
telescopes. The reason is that telescopes are sensitive to showers of energies up to four orders 
of magnitude below the average energy of cosmogenic neutrinos. The excess of hadronic 
showers at IceCube could be observed for M D up to 5 TeV. 
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